SLAC-PUB-7634 
UPR-788-T 



Bounds on Supersymmetry from 
Electroweak Precision Analysis 



Jens Erler* 

Department of Physics and Astronomy 
University of Pennsylvania 
Philadelphia, PA 19104 

Damien M. Pierce^ 

Stanford Linear Accelerator Center 
Stanford University 
Stanford, CA 94309 



Abstract 

The Standard Model global fit to precision data is excellent. The Minimal Supersym- 
metric Standard Model can also fit the data well, though not as well as the Standard 
Model. At best, supersymmetric contributions either decouple or only slightly decrease 
the total x 2 ] a t the expense of decreasing the number of degrees of freedom. In general, 
regions of parameter space with large supersymmetric corrections from light superpart- 
ners are associated with poor fits to the data. We contrast results of a simple (oblique) 
approximation with full one-loop results, and show that for the most important observ- 
ables the non-oblique corrections can be larger than the oblique corrections, and must 
be taken into account. We elucidate the regions of parameter space in both gravity- and 
gauge-mediated models which are excluded. Significant regions of parameter space are 
excluded, especially with positive supersymmetric mass parameter /i. We give a com- 
plete listing of the bounds on all the superpartner and Higgs boson masses. For either 
sign of fi, and for all supersymmetric models considered, we set a lower limit on the 
mass of the lightest CP-even Higgs scalar, mn > 78 GeV. Also, the first and second 
generation squark masses are constrained to be above 280 (325) GeV in the supergravity 
(gauge-mediated) model. 



* e-mail: erler@langacker.hep.upenn.edu 
t e-mail: pierce@slac.stanford.edu 



1 Introduction 



For more than a decade, ideas involving supersymmetry (SUSY) have been among the most 
popular extensions of the Standard Model (SM), having the greatest potential to solve its 
shortcomings such as the gauge hierarchy problem and the lack of a quantized version of grav- 
ity Indeed, all superstring theories necessarily contain quantum gravity and supersymmetry. 
Moreover, very recently, using arguments based on various dualities and supersymmetry, all 
superstring theories as well as 11 dimensional supergravity seem to be connected nonpertur- 
batively |1[, suggesting a unified theory in which supersymmetry is one of the indispensable 
key elements. The idea of super symmetric unification is further supported by the observation 
of gauge coupling unification [0 two orders of magnitude below the reduced Planck scale, 
Mp/-y/87r ~ 2 x 10 18 GeV, the natural supergravity scale. 

Taken together this implies a strong motivation to investigate the phenomenological conse- 
quences of low energy supersymmetry. In this paper we present a systematic study of precision 
observables in the minimal supersymmetric standard model (MSSM) in an attempt to find 
favorable and/or excluded regions in SUSY parameter space. We do this in two scenarios for 
how supersymmetry breaking is conveyed to the observable sector. 

In the "minimal supergravity" model || supersymmetry breaking is transmitted from 
a hidden sector to the observable sector via gravitational interactions. Supersymmetry is 
spontaneously broken in the hidden sector at a scale A H. In popular models which invoke 
gaugino condensation ||, A ~ (MpMsusY /871-) 1 / 3 ~ 10 13 GeV. Since gravity is flavor blind, it 
is assumed that the explicit soft breaking terms in the observable sector are universal at the 
supergravity scale. 

In simple models of gauge mediation || there is a supersymmetry breaking sector which 
gives rise to both an F-term, Fx, and a vacuum expectation value, X, of a standard model 
singlet field. This field is coupled to the vector-like "messenger fields" , M, through a superpo- 
tential interaction of the form XX MM. The generated soft mass term for a given superpartner 
is proportional to A = F x /X ~ 4nM z /a ~ 10 5 GeV, and grows with the square of its gauge 
couplings. Therefore sleptons are much lighter than squarks and the gluino is heavier than 
charginos and neutralinos. 

In both cases the generated soft terms are sufficiently flavor universal that all flavor chang- 
ing neutral current (FCNC) effects are adequately suppressed. However, in the minimal su- 
pergravity model it must be implicitly assumed that universality is not spoiled by the Kahler 
potential. 

Recently, a number of articles examined whether supersymmetry has the potential to 
describe the data better than the SM. The conclusions in the affirmative were mainly driven 
by Rj, (see section Q), which at times was more than 3cr higher than SM expectations. Loops 
involving light charginos and top squarks could account for a large part of the discrepancy 
for low tan (3, provided the effect was not canceled by charged-Higgs loops. For large tan /3, 
large shifts could be obtained from loops containing bottom quarks and light neutral Higgs 
bosons, and to a lesser extent from neutralino/bottom-squark loops. However, more recent 
analyses find the measured value of closer to the SM prediction (the discrepancy is 1.3<r), 
and at the same time direct limits on superpartner masses have increased. Now, even if one is 
able to find a region of parameter space where the Rb discrepancy is alleviated, the decrease 
in the overall \ 2 is 110 1 significant. Also, in those regions of parameter space one typically 
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finds that the discrepancy in A FB (b) is made worse. In the specific high-scale models of 
supersymmetry breaking we consider in this paper, the largest possible shift in Rb is less than 
la. It follows that Rb now plays a much smaller role, and supersymmetric models can no 
longer yield significantly smaller values of y 2 than the SM. 

Therefore, in this work we take a different point of viewQ and focus on elucidating the 
excluded regions of supersymmetry parameter space. We present a complete one- loop analysis 
of supersymmetry, combined with a state of the art SM calculation. Input data are as of August 
1997. 

In Section |2| we give a short overview of the inputs and observables we use. For reference we 
perform a global fit to the SM, and comment on some of the deviating observables. Section 
describes the parameter spaces of the supersymmetric models in more detail, and reviews 
recent direct limits on superpartner masses. We compare the oblique approximation \ V2\ with 
the full calculation in some detail. We present our results in Section £| and our conclusions in 
Section [S]. 



2 Overview 

Within the SM we perform a global fit to a total of 31 observables listed below. We include 



full one-loop radiative corrections full QCD corrections up to 0(al); higher order QCD 
corrections when enhanced by beta function effects; mixed electroweak/QCD corrections of 
0(aa s ) with the exception of non-leading special vertex corrections to Z — > bb decaysQ; and 
0(a 2 ) corrections when enhanced, e.g. by large top mass effects or large logarithms. 

The Fermi constant, G M , and the fine structure constant, a, are taken as fixed inputs. 
The five fit parameters are the strong coupling constant, a s ; the contribution of the five light 
quarks to the photon vacuum polarization function, Aa^; and the masses of the Z-boson, 
Mz, the Higgs scalar, M#, and the top quark, m t . Alternatively, one may fix Mz as well, since 
its relative error is now comparable to that of G^. We have chosen to leave it free, because the 
Z-mass measurement is correlated with other observables. In practice, the difference between 
the two treatments is numerically insignificant. We have organized the measurements into 
seven groups: 

• 9 lineshape observables 

The mutually correlated LEP observables [16[ are determined from a common fit to the 



Z lineshape and the leptonic forward-backward asymmetries. They include the Z-boson 
pole mass Mz, the total Z-width Tz, the hadronic peak cross section, 

_ 127rr 

and, for each lepton flavor, £ = e, //, r, the ratios Ri = Fh a d/r«, and the pole asymme- 



1 First results of this kind of analysis were presented in Ref. ||. 

2 These have been calculated very recently and shown to be quite small compared to the analogous correc- 
tions to the lighter quark vertices |15[ . 
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tries, A FB (£). T x denotes the partial Z width into x. Defining 

1 - 4Q, sin 2 9 f cS 



A 



f 



1-4Q /S in 2 e f eS + 8Qj sin 4 0, 
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where sin 2 9( s is the effective weak mixing angle for fermion / at the Z scale, we have 



A*»(J) 



;A e A f . 



A T , its forward-backward asymmetry, V (r 



3 further LEP asymmetries 

These are the r polarization, V{j) 
A e , and the hadronic charge asymmetry, (Q FB ), which is quoted as a measurement of 

6 heavy flavor observables 

The mutually correlated heavy flavor observables from LEP and SLC [[UJ are the ratios 
R q = rgg/Thad; the forward-backward pole asymmetries, A FB (q) (LEP); and the com- 
bined left-right forward-backward asymmetries, A F ^(q) = A q (SLC), each for q — b,c. 

3 further SLD asymmetries 

These are the very precise measurement of the left-right asymmetry for hadronic final 
states []T7 ] , A^had) = A e ; the analogous Alr(h,t) = A e , which is obtained from a 
common fit including polarized Bhabha scattering |Tj|; the left-right forward-backward 
asymmetries [IS|, A FB (£) = A^ for I = \x or r; and the hadronic charge flow asym- 
metry |19 |, Aq = A e , which is basically the ratio of weighted forward-backward and 
left-right forward-backward asymmetries^. We combine the 3 values of A e from SLD 
into one number, 



A P = 0.1548 ±0.0033 



2 pole masses 



We combine the VT-mass measurements from CDF [20], D0 and UA2 [22|, M\y(pp) 



80.409 ± 0.090 GeV, with the one from LEP 2 @, M W (LEP) = 80.481 ± 0.140 GeV, 
obtaining 

Mw/(world) = 80.430 ± 0.076 GeV. 



Combining the results of all top decay channels at CDF |23) and D0 

m t = 175 ± 5 GeV. 

We fix the scale invariant MS masses for the heavy quarks, mb(jn b ) 
m r (m r ) = 1.30 GeV. 



we find 



4.33 GeV and 



3 Sincc Aifl(had) is proportional to the electron beam polarization, V e , while Aq is inversely proportional, 
the geometric mean of the two determinations yields A e independently of V e - The result jl9| is A e = 0.1574 ± 
0.0208, or equivalents, sin 2 6 e cS = 0.2302 ± 0.0027. 
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6 low energy observables 

The weak charge from atomic parity violation (APV) in Tl has been measured by groups 
in Oxford |5[ and Seattle 0, Q W ( 205 T\) = -114.77 ± 1.23 ± 3.44, and in Cs by the 
Boulder group [p7fl , Qw( 133 Cs) = —72.11 ± 0.27 ± 0.89, where the first errors are ex- 
perimental and the second theoretical. The recent result of the deep inelastic scattering 
(DIS) experiment of CCFR |28| is combined with the measurements of CDHS |29] and 



CHARM [pcfl , yielding k = 0.5805 ± 0.0039, where k is a linear combination of effective 
4-Fermi operator coefficients. Similarly, v^e scattering experiments, dominated by the 
recent CHARM II results ||31|| , yield determinations of leptonic 4-Fermi operator coef- 



ficients, ffi = -0.041 ± 0.015 and g v £ = -0.507 ± 0.014. The CLEO measurement of 
B — > X s j |33[] yields the 90% confidence interval, 



1 x 10~ 4 < B(B -> X sl ) < 4.2 x 10" 4 . 



When a variable is defined in a finite interval (here between and 1), and a bound- 
ary of the interval is within a few standard deviations from the central value, a vari- 
able transformation should be performed such that the new variable is defined on the 
whole real axis. Arguments analogous to the ones supporting a flat prior distribution 
in log Mh rather than in Mh (see next section) lead to the logistic transformation?], 
lg(x) = ln(x/(l — x)) ||34|1 . This results in a more Gaussian shaped distribution, and in 
our case, in a more (less) conservative assessment of the upper (lower) error bar. The 
b — > S7 amplitude can receive very large contributions from the leading-order supersym- 
metric loops, and the next-to-leading order corrections (which have not been calculated) 
could be important. The uncertainty in the large and positive MSSM contributions is 
taken into account by the more conservative treatment of the upper error bar. Addi- 
tionally, because of these uncertainties, we have chosen the more conservative of the two 
error estimates in Ref. |33fl , which is obtained by adding the statistical and systematic 



errors linearly. Compared to this experimental error, the theoretical error due to QCD 
uncertainties in the SM prediction |35| can be neglected. 



2 gauge couplings 

We use the constraint AaJ^j = 0.02817 ± 0.00062 @. We also include the external 



constraint a s = 0.118 ± 0.003, which we obtained by combining non-Z lineshape deter- 
minations of a s |57|. In this case we have doubled all theoretical errors to account for 
ignored correlations, and because these are by far the most difficult to estimate and very 
non- Gaussian in nature. 

The best fit values of the SM input parameters are shown in Table |l| The x 2 /d.o.f. of the 
fit is 26.6/26, corresponding to an almost perfect fit with a goodness (probability of larger \ 2 ) 
of 43%. Our results differ from the ones quoted by the LEP Electroweak Working Group |16| . 



The differences can be traced to a different treatment of radiative corrections and to a slightly 
different and more recent data set. Our results are in close correspondence with the very recent 
update of the Particle Data Group [|T3J; one major difference is the external a s constraint. 



4 Equivalcntly, one may continue to use the old variable and include a nontrivial weight factor (the Jacobian) 
l/[x(l-x)]. 
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m t 


172 ± 5 GeV 


M H 


66t|| GeV 


M z 


91.1867 ±0.0020 GeV 


a s [M z ) 


0.1195 ±0.0022 




0.02817 ±0.00064 



Table 1: Results for the SM fit parameters. 



Also, here M# is allowed as a free parameter, in which case the SM best fit value of Mh is 
below the lower limit from direct searches. In the MSSM the bound is relaxed to nih > 60 
GeV [38] , consistent with the central value in Table [l|. 



We list the measured values of all the observables in Table |2], together with their global 
best fit values. The pull, defined as the difference between the measurement and the fit result 
divided by the experimental error, is also shown. The agreement is excellent. The two largest 
discrepancies are only at the 2a level. 

Based on all data from 1992-1996, the left-right asymmetry, Ai^(had) = 0.1550±0.0034, is 
now closer to the Standard Model expectation of 0.147 ±0.002 than previously |17]]. However, 
because of the smaller error, is still significantly above the Standard Model prediction. 
There is also an experimental difference of ~ 1.9cr between the SLD value of Ag (SLD) = 
0.1547 ± 0.0032 from all A LR and Af| (£) data on one hand, and the LEP value A e (LEP) = 
0.1461 ± 0.0033 obtained from A FB (£), A e (V T ), and A T (V T ) on the other hand, in both cases 
assuming lepton-family universality. 

Relaxing universality, one can extract the various Af in a model independent way. First, 
we combine A e from A FB (e) and V fb (t) with the value from the SLC to obtain A e = 0.1518± 
0.0029. This value for A e can now be used to extract the other Af from forward-backward 
asymmetries. Looking at Table one sees that all 3 observables which are sensitive to the r- 
vertex, namely A fb (t), V(t), and A f ^(t), deviate by la or more. However, when combined, 
one finds A T = 0.1448 ±0.0059, in excellent agreement with the SM value of 0.147 ± 0.002. At 
present, it seems more likely that we are facing an experimental discrepancy. Using A FB (f) 
from LEP and A FB (f) from the SLC for f = fi, b, and c, we find that A^ and A c are consistent 
with each other and the SM. On the other hand, while the Aj, determinations are also consistent 
with each other, their combined value, Af, = 0.872 ± 0.024, deviates 2.6<r from the SM value 
of 0.935. A similar analysis assuming lepton universality yields |T6| A b = 0.877 ± 0.023 which 
is 2.5cr low. One can also obtain A e = 0.148 ± 0.002 from a global fit excluding the quark 
sector^, and derive Ab = 0.888 ± 0.022, still a 2.1a deviation. We can therefore conclude 
that the data show an anomaly in A&, independently of exactly how it is extracted. However, 
this deviation of about 6% cannot be due to supersymmetric loops since a 30% correction to 
k b (defined by sin 2 6 b cS = k b sin 2 9^) would be necessary to account for the central value of 
A b . Only a new contribution at tree level which does not contradict Rb (including the off- 
peak Rb measurements by DELPHI |39|]), can conceivably account for such a low Ab 4(|. In 



particular, a r-sneutrino with mass close to M z and R-parity violating couplings is a natural 
candidate ||41]. Of course, supersymmetry can also affect A FB (b) through the weak mixing 



'Here we have fixed Mh = Mz because otherwise the Higgs mass is driven to very low (excluded) values. 
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angle in the A e factor. Such a correction can not, however, help to resolve the LEP/SLC 
discrepancy. 

This discussion and the remarks about Rb in the introduction show that we cannot expect 
that the inclusion of supersymmetric loops will significantly improve the goodness of the fit 
to electroweak data. On the contrary, there are regions in parameter space where the fit is 
quite poor, and those regions must be ruled out. Here we do not test the supersymmetry 
hypothesis. Rather, the primary goal of our analysis must be to identify the excluded regions 
of supersymmetry parameter space. 



3 Supersymmetry 

In gravity-mediated models supersymmetry is broken in a hidden sector and the breaking 
is communicated to the observable sector through supergravity effects. This transmission is 
suppressed by inverse powers of the Planck mass. The assumption of a flat Kahler metric for 
the fields in the observable sector (an assumption which can be called into question) yields 
a universal scalar mass, M , a universal gaugino mass, M1/2, and a universal trilinear scalar 
coupling, Aq, at the Planck scale |ji2| . In the minimal supergravity model we impose this 
universality at the scale where the (U(l) and SU(2)) gauge couplings unify, Mgut ~2x 10 16 
GeV. The running of the soft parameters from Mp to Mqut should in general be taken into 



account [53], but it is model dependent, so it is ignored in the minimal supergravity model. 

From the inputs M Z} a, Aa^ d , and a S} and the fermion masses m t , m b} and m T , we 
determine the full one-loop values of the gauge and Yukawa couplings at the Z-scale. We then 
solve the set of coupled two-loop renormalization group equations (RGE's) for the couplings 
and soft parameters by iteration. This set of differential equations has two-sided boundary 
conditions: the hard parameters (gauge and Yukawa couplings) are fixed at the Z-scale and 
the soft parameters are fixed at the unification scale. In each iteration we impose electroweak 
symmetry breaking at the squark scale. This occurs naturally, as the large top-quark Yukawa 
coupling drives the Higgs mass-squared m 2 ^ negative. The electroweak symmetry breaking 
conditions allow us to solve for the Higgsino mass-squared, /x 2 , and the CP-odd Higgs boson 
mass-squared, m\. These are determined as functions of tan /3, the ratio of vacuum expectation 
values, to full one- loop order ||10|| . To summarize, the minimal supergravity model parameter 
space is specified by 

M , M1/2, A , tan/?, and sgn(/i) . (1) 

Just as in minimal supergravity, simple models of gauge-mediated supersymmetry breaking 
also introduce a hidden sector in which supersymmetry is dynamically broken, so that for one 
(or more) hidden sector fields X, Fx 7^ 0. Here, however, the messenger fields which couple to 
X are assumed to be charged under the SM gauge group. As a result, supersymmetry break- 
ing is transmitted through messenger loop effects, and observable fields receive soft masses 
roughly proportional to the square of their gauge couplings. It follows that scalars with iden- 
tical quantum numbers are mass degenerate, and FCNC effects are sufficiently suppressed as 
long as the messenger scale M = XX < 10 15 GeV [44]. The key parameter here is A = Fx/X, 



as this parameter sets the scale of the entire superpartner spectrum. The messenger scale 
M determines the amount of RGE evolution, and therefore enters only logarithmically. As 
before, radiative electroweak symmetry breaking occurs rather generically over the parameter 
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space. It is generally difficult to arrange for the right magnitude of fi and the soft supersym- 



metry breaking parameter B, but solutions to this problem exist ||45|| . We assume that the 
mechanism which is responsible for generating B and [i does not at the same time give extra 
contributions to the soft squark and slepton masses. The superpartner spectrum depends on 
the representation(s) of the messenger sector. We consider the cases where the messenger 
sector is comprised of a 5 + 5 pair of fields, or a 10 + 10 pair[]. For the purpose of our analysis 
a 10 + 10 messenger sector is equivalent to a model with three 5 + 5 pairs. In the minimal 
gauge-mediated model the supersymmetry parameter space includes 

M, A, tan/3, and sgn(/i), (2) 

with A < M to ensure the absence of gauge symmetry breaking in the messenger sector. 

In either supersymmetry model, we proceed by randomly picking a point in its parameter 
space, Eq. ([!]) or Eq. @. For an initial set of values for the standard parameters , Mz, m t , a s 

(5) 

and Aajjad, we then solve the model, i.e. we self-consistently evaluate the soft parameters at 
the weak scale, and compute the physical superpartner and Higgs-boson masses, couplings and 



mixings. This enables us to compute the full supersymmetric one-loop corrections |I0], [46|, |47 
to each observable^], and to perform a y 2 minimization with respect to the standard parameters. 
We iterate the solution of the supersymmetric model for the best fit values of the standard 
parameters until this process has converged. On the last iteration we require that the full 
one-loop electroweak symmetry breaking conditions are satisfied, and impose that the Yukawa 
couplings remain perturbative (A 2 /47r < 1) up to M G tjt- 

Finally, we apply the direct lower limits on superpartner and Higgs boson masses. Searches 
at LEP constrain the mass of the CP-odd Higgs boson, A, to be larger than 62.5 GeV, and 
that of the lightest CP-even Higgs scalar, h, to be heavier than 62.5 GeV for small m^, or 70.7 
GeV for large |j48| . There are mass limits from LEP 2 on charginos {m^+ > 91 GeV) 



selectrons (m g > 76 GeV), smuons (m^ > 59 GeV), staus (m f > 53 GeV) pD[ , and stops 
{nil > 65 GeV) |5l| provided these superpartners are not nearly degenerate with the lightest 
neutralino, Xv The chargino limit is slightly degraded if the electron sneutrino is light. The 
top squark limit is for the case when it decouples from the Z, and increases to 72 GeV for 
no left-right mixing. For a relatively light Xi, D0 excludes top squarks up to 95 GeV [52|. 



In the gauge-mediated models the next-to-lightest supersymmetric particle (NLSP) can be a 
slepton. We assume in these models that the NLSP decays outside the detector, in which case 
the search limits for stable charged particles apply, which are > 67 (69) GeV for right- 
handed (left-handed) sleptons f53 |. The for us relevant CDF searches |54| suggest for first 



and second generation squarks, nig > 210 GeV, and for gluinos, rrig > 173, but these limits 
do not apply uniformly over parameter space. Therefore we impose relaxed limits of 200 and 
160 GeV, respectively. However, it turns out that in the models we consider the direct limits 
on squark and gluino masses are for the most part irrelevant. Once the other mass limits are 
taken into account, the smallest possible squark (gluino) masses vary from 210 (230) to 370 
(330) GeV, depending on the model. 

6 The messenger fields are assumed to come in complete SU(5) representations to preserve the successful 
gauge coupling unification. An actual SU (5) symmetry, however, is not assumed. 

7 For the low-energy observables from atomic parity violation, deep inelastic scattering, and neutrino scatter- 
ing, we work in the oblique approximation. For these observables the non-oblique supersymmetric corrections 
are expected to be negligible relative to the experimental errors. 
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At each point in our scan over parameter space we apply the constraints from electroweak 
symmetry breaking, Yukawa perturbativity, and the direct searches. If any of these checks 
fail, the point is disregarded. We refer to a point in the supersymmetry parameter space 
which passes these tests as a quaint point. The set of points which pass comprise the quaint 
parameter space. 

We now discuss the interpretation of the resulting y 2 value for a given quaint point. It 
is important to notice that the overall x 2 value is relevant only for hypothesis testing. As 
noted before, this is not our goal. We derive our confidence that the study of low-energy 
supersymmetry is worthwhile from other (more indirect) arguments. Rather, we will address 
the question of preferred or disfavored quaint regions. Before we continue to defend our 
procedure, we would like to compare our situation with the SM Higgs mass. 

In standard analyses, the quoted limits on Mh are independent of the \ 2 value of the SM 
best fit. A standard (Bayesian) procedure is to determine the posterior probability distribution 
function (posterior p.d.f.) as the product of the prior p.d.f. (prior) and the likelihood function 
(here ~ exp(— % 2 /2)). The choice of a prior is always made, but is oftentimes obscured. In 
the case of Mh, a non-informative prior is desired, or perhaps an informative prior which 
accounts for the direct lower limit, and is non-informative otherwise. The choice of a suitable 
prior is not unique^. A non- informative prior is by no means necessarily flat (constant); in 
general a flat prior ceases to be constant after a variable transformation introduces a non-trivial 
Jacobian. Clearly, such a variable transformation cannot change the results of a statistical 
analysis unless the Jacobian is dropped. In the case of Mh, one often chooses a prior which is 



flat in log Mh. There are various ways to justify this choice |34|| . One rationale is that a flat 
distribution is most natural for a variable defined over the real numbers. This is the case for 
log Mh but not M H . Another complication arises since flat priors in both logM^- and M H are 
improper, i.e. their integrals are divergent and cannot be normalized. This is not necessarily 
problematic as long as the likelihood function is sufficiently convergent so that the posterior 
p.d.f. is normalizable. For Mh this is the case, due to the non- decoupling property of the SM 
Higgs boson corrections. 

In the case of supersymmetry, however, the posterior p.d.f. remains improper, because 
superpartners and the heavy Higgs doublet decouple, and the likelihood function is asymptot- 
ically finite. To cope with this problem one might 

• cut superparticle masses or soft parameters off somewhere in the TeV region, arguing 
that higher scales are unnatural and incompatible with the primary motivations for 
low-energy supersymmetry such as the hierarchy problem. In other words one would 
introduce an informative prior. However, conclusions would directly depend on the 
precise choice of the cut-off scale, i.e. the prior; 

• choose a more sophisticated informative prior which utilizes fine-tuning criteria in order 
to produce a proper posterior p.d.f. While aesthetically appealing in particular from 



a Bayesian viewpoint ||55||, the implementation would be rather involved and computa- 
tionally expensive. Also the establishment of an appropriate fine-tuning criteria is not 
unique [EGA; 



8 The sensitivity of the posterior p.d.f. to the (non-informative) prior diminishes rapidly with the inclusion 
of more data. 
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• exclude all quaint points which have a Ax 2 = x 2 ~ Xmin > 3.84, where Xmin i s the 
minimum value of x 2 m the supersymmetry model under consideration. In a univariate 
Gaussian situation this procedure corresponds to the 95% CL. It lacks the conceptual 
rigor of Bayesian statistics, but has the advantage of being independent of the specifica- 
tion of prior information. We will use this criterion in our analysis. 

Before showing the results of the full one-loop analysis, we discuss the results obtained in 
the oblique approximation. For some observables the dominant supersymmetric corrections 
are of oblique type, i.e. they arise from vector-boson self-energy diagrams and can be described 
by the S, T and U parameters^ fl2| . These parameters can only be extracted from the data 
when Mh is fixed to a reference value. The following approximate formulas reproduce S, T 
and U within ±0.01 for 60 GeV < M H < 1000 GeV: 

S = -0.12 ± 0.14 + 0.04(i^ _ 1} _ 13 logiQ _a^_ ; 

T = -0.10 ± 0.14 + 0.08(i^ - 1) + 0.31 log 10 (3) 

U = +0.20 ± 0.24 - 0.02( 10 °^ cV - 1). 

We see that all 3 oblique parameters are consistent with zero within la. Now consider the 
supersymmetric contributions to 5*, T and U. Supersymmetry contributes to the S parameter 
with either sign, but that changes the overall x 2 only by relatively small amounts typically 
of 0(1). Contributions to T and U are positive definite, but U is increased by less than 0.1, 
which is much less than the present uncertainty. On the other hand, contributions to T can 
be significant, and since positive, Eq. (|3|) shows that strong constraints can be expected to be 
associated with T. 

The general oblique fit in Eqs. (|3|) is not very suitable for supersymmetry. It is more 
appropriate to impose the constraints T > 0, U > 0, and 60 GeV < Mh < 150 GeV, and to 
deviate from the oblique approximation in allowing an extra parameter, 75, defined through [5T5] 



T(Z — > 66) = T°(Z — > 66) (1 + 7&), in order to account for the possibility of supersymmetric 
Zbb vertex corrections. The result is^ 

s = -o.o9±g:ig, 

T < +0.06, 

u = +o.irt° : f 7 , (4) 

7 6 = 0.006 ± 0.004. 

If we redo this fit with % = 0, the strong constraint, T < 0.06 is relaxed to T < 0.08. 
Eq. ([|) shows that the T parameter is the most important of the three oblique parameters 
in determining the goodness of the fit. We show the various supersymmetric contributions to 
the T parameter in Fig. [I] in the supergravity model. The shaded regions are determined by 
randomly choosing O(10 5 ) points in parameter space ([!]) to find 50,000 quaint points. We see 
that the maximal contributions from each of the three superpartner sectors, gauginos, third 
generation squarks, and sleptons are significant (£T max ~ 0.06, 0.12, and 0.08, respectively). 
The Higgs sector (and the first two generations of squarks) contributes at a lower level (<5T max ~ 
0.015). This pattern holds for the S and U parameters as well. 

9 We use the definitions from reference JoTj , but redefine them to include contributions from new physics 
only|8). 

10 This fit does not include the constraint from B(B — » Xg^). 
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Figure 1: Supersymmetric contributions to the T parameter in the minimal supergravity 
model. Shown are (a) the chargino/neutralino contributions vs. the light chargino mass; (b) 
the stop/sbottom contributions vs. the heavy stop mass; (c) the slepton contributions vs. the 
left-handed selectron mass; and (d) the Higgs sector contributions vs. the CP-odd Higgs boson 
mass. Masses are in GeV. 



We show the total contribution from all superpartners and extra Higgs bosons to S, T 
and U in Fig. 0. Here we see that these contributions are approximately in the range (—0.05 
to +0.08), (0 to 0.2), and (0 to 0.09), respectively, for S, T and U, and the decoupling 
of supersymmetry from these parameters is evident. These contributions are significantly 
limited by the recent, more stringent collider bounds on new particle masses discussed above. 
In the gauge-mediated models, the maximum contribution to the oblique parameters is even 
further reduced. We illustrate this in Fig. |3|, which shows the allowed regions in the S-T plane 
for the three models. We also show the best fit values of S and T, and the 68% and 95% 
contours, determined with Mh fixed at 100 GeV. It is important to realize that the difference 
in x 2 between the SM minimum and the minimum allowing freely varying oblique parameters 
is only one unit. Hence, even with the best possible addition of oblique corrections, the \ 2 
will not improve significantly with respect to the SM. Moreover, Fig. ^| nicely illustrates two 
points. First, the corrections from supersymmetry are quite limited in magnitude, and second, 
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Figure 2: Supersymmetric contributions to the S, T, and U parameters vs. the light chargino 
mass, in the supergravity model. 



they tend to worsen the fit compared to the SM. 

Judging from Fig. 3, it appears that no or very little of the supersymmetric parameter 
space can be excluded. However, the oblique approximation is just that, an approximation. 
As we will show, it cannot be used to reliably indicate whether parameter space is excluded. 
The oblique approximation neglects higher order terms in the derivative expansion, as well as 
the non-universal corrections: the vertex corrections, wave-function renormalization, and box 
diagrams. We illustrate the relative magnitude of the oblique and non-oblique corrections to 
many of the observables in Fig. |4], in the supergravity model. Here, we show the full range 
of the shifts of each observable due to the shift in each oblique parameter]^], labeled S, T, 
U, and their sumPI labeled obi. For comparison, just below the range of the shift from the 
total oblique correction we indicate the non-oblique corrections at the endpoints of that range. 
Therefore one can compare the oblique and non-oblique corrections at the points of maximal 
positive and negative oblique shift. We also show the range of the non-oblique corrections 
in our scan, labeled n-obl, and just below that, for comparison, the corresponding oblique 
correction at the points of maximal positive and negative non-oblique shift. Finally, we show 
the range of the full one-loop corrections, labeled full, and at the points of maximal positive 
and negative full correction we indicate the relative sizes of the oblique and non-oblique parts. 
All these shifts are shown in units of standard deviation, i.e. we divide each shift by the 
experimental error. The APV observables, the neutrino scattering observables, as well as 
A^|f(/i), A^ft ip), A£jf(c), and R c , all receive very small shifts, and are not shown. We show 
the same plots for the 5 + 5 and 10 + 10 gauge mediation models in Figs. |^ and ^. From 
these figures we see that, in all three models, the non-oblique corrections can be larger than 
the oblique corrections for every observable except My/ (and Y z in the supergravity model). 

n Only Mw receives non-negligible contributions from the U parameter. In particular, the DIS variable k 
is insensitive to U, although DIS results are frequently quoted in terms of sin 2 9w = 1 — Myy/M%. Clearly, in 
the presence of new physics such quotes can be misleading. 

12 The shift from the maximal total oblique correction is not necessarily the sum of the shifts due to the 
individual maximal oblique corrections, because the individual corrections need not find their maximal values 
at the same point in parameter space. 
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Figure 3: Supersymmetric contributions to the S and T parameters. Shown are the results for 
(a) the minimal supergravity model; (b) the 5 + 5 gauge-mediated model; and (c) the 10 + 10 
gauge-mediated model. The best fit value (+) and 68 and 95% CL contours of S and T are 
also shown, for M H = 100 GeV. Note that by definition S = T = in the SM. 



Hence it is essential to include the full corrections in this type of analysis. 

Figs. |]-|6] also show that the largest corrections occur in the supergravity model, and the 
smallest in the 5 + 5 gauge-mediated model. In each model the gaugino masses obey the 
usual GUT relation, i.e. oc c^. For a given gaugino mass scale, the scalar spectra (and 
the Higgsino mass) in the three models vary. The model with the lightest scalars (for a fixed 
gaugino mass) is the supergravity model, in the region M <C Mi/ 2 . It is in this region (actually 
M ~ Mi/2/2) that the largest corrections occur. The 5 + 5 gauge-mediated model has the 
heaviest scalar spectrum, and hence the smallest corrections (Fig. ||). The scalar masses can 
be made relatively lighter in the gauge-mediated case if the effective number of 5 + 5 pairs, 
n| , is increased. The 10 + 10 model we consider corresponds to n e ^ = 3, and hence the 
corrections are larger in the 10 + 10 case (Fig. ||). However, they are still not as large as in 
the supergravity model. We expect the magnitude of the corrections in the case0 n| = 4 
to match quite closely with the supergravity results, as in this case the two spectra can have 
significant overlap (see Figure 10, Table 1, and the corresponding discussion in Ref. [B[J, for 
an example). Figs. [5] and ^ show that insensitive observables in the supergravity model are 
likewise insensitive in the 5 + 5 and 10 + 10 gauge mediation models. 

Figure [7] shows, for each model, the range of the best fit predictions of each observable 
found in our scan over quaint parameter space, in units of standard deviation. We also indicate 
the best fit value of each observable in the SM. Note that the range of each observable in 
the supersymmetric models always includes the SM best fit value, since the supersymmetric 
corrections decouple, and the range of the supersymmetric light Higgs boson mass includes 
the best fit SM Higgs mass. The best fit values of the observables which are most sensitive to 
supersymmetric contributions, Tz, A FB (b), and the SLC measurement of A e , can vary over a 
range of ~ la. This range for the best fit predictions is smaller than the actual contributions 
from supersymmetry shown in Figs. which can be larger than 2a. This is because the 

13 This corresponds to an SO(10) model with a 16 + 16 pair of fields. 
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standard parameters adjust themselves to minimize the overall x 2 . The sensitivity of Tz, 
A FB (b), and A e are due to the very high precision with which these observables are measured. 
Shifts in My/ are somewhat smaller, but still important. Taken together, these large shifts do 
not improve but rather worsen the fit. They give rise to regions in parameter space which can 
be excluded. We describe these excluded regions in the next section. Not shown in Fig. |7| are 
correlations among the supersymmetric corrections. These correlations figure importantly in 
the fit. For example, in the region of parameter space where the supersymmetric prediction 
for R b is closest to the measured value, the discrepancy in A FB {b) is much worse. 



4 Results 

The allowed and excluded quaint regions in a variety of two-dimensional parameter subspaces 
are shown in Figs. [§] and |9| for fi < 0, and in Figs. [II], [IT], and |T2] for fi > 0. These plots are the 



main result of our analysis. Each row of these figures shows the results for the supergravity 
model, the 5 + 5 gauge-mediated model, and the 10 + 10 gauge-mediated model, respectively 
from left to right. We show with a solid line the region of parameter space where it is 
possible to find A% 2 < 3.84. With a dashed line we indicate the region where it is possible 
to find Ax 2 > 3.84. Hence, those regions enclosed by a dashed line but not by a solid line 
are excluded, independently of the value of any other parameter. Of course, in the regions 
bounded by a dashed line and a solid line there are both allowed and excluded points. In 
some cases we plot the shift 5 from supersymmetry rather than the value of the observable. 
We show more results in the /i > case because of the larger excluded region, as discussed 
below. In the following discussion, we use the abbreviations SUGRA, GM 5 , and GM 10 to refer 
to the minimal supergravity, 5 + 5 and 10 + 10 gauge-mediated models. Also, GM refers to 
both the 5 + 5 and 10 + 10 models, and limits listed as three consecutive numbers refer to the 
SUGRA, GM 5 , and GMi models, respectively. 

Figs. H (a-c) and [10] (j-1) demonstrate the importance of the B(B —>■ A s 7) observable for 
our results^]. When included much larger values of A% 2 are possible. Large parts of parameter 
space are excluded due to unacceptably large positive or negative contributions to the b — ► sj 
amplitude, especially with /i > 0. The most important corrections are due to chargino/top- 
squark loops. These contributions are proportional to /itan/3, at large tan (3. Hence, for 
\i > 0, they constructively add to the SM and charged Higgs amplitudes, and very large total 
amplitudes can result over a wide region of parameter space. For large negative /z, the chargino 
amplitude destructively interferes with the SM and charged Higgs amplitudes, and the full 
one-loop amplitude can vanish. However, this cancellation occurs only in a relatively small 
region of parameter space. Either case can result in a very large contribution to A% 2 . Because 
of the strong impact of this observable we have chosen a more conservative error assignment 
for it (see Section |2[). 

As a first example of how the /i > region of parameter space is more severely constrained 
than the fi < region, consider the input parameters in the supergravity model. In the fi > 

14 Since the next-to-leading order calculation of B(B — > X s j) is not available in the MSSM, we use an 
improved leading order formula in the Standard Model calculation and include the leading order super- 
symmetric corrections p2[. 
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case, the no-scale model |k| (i.e. M = 0) is seen to be excluded^ Specifically, we find the 
independent limits M > 9 GeV and M 1/2 > 105 GeV (Fig. 0(a)). In the /i < case, by 
contrast, no limit on M can be set, and the limit, M l / 2 Z 120 GeV, is not improved by our 
analysis. 

We should stress at this point that, for either sign of fi, there is some parameter space 
excluded even at the very largest values of M and Mw 2 that we consider (1 and 0.5 TeV, 
respectively) . Here we emphasize the bounds that are valid independently of the values of any 
other parameters. This is what we mean by "excluded regions". For emphasis we sometimes 
refer to these regions as "absolutely excluded" . 

The input parameter which sets the scale of the superpartner spectrum in the gauge- 
mediated models, A, is constrained in the \x > case to be larger than 36 (14) TeV in the 
GM 5 (GMio) model, whereas the entire (quaint) parameter space extends down to 28 (10) 
TeV (Figs. [10] (b-c)). Again, the corresponding bound with fi < of 41 (13) TeV could not 
be improved significantly. 

Similarly, Figs. [□] (d-f) show that for \i > 0, tan/5 < 36, 45 and 39, respectively, for the 
three models. In contrast, for fi < values of tan/? as large as 65 are still allowed. 

Figs. [K] (d-f) illustrate how the oblique parameters can serve as useful barometers in 
parametrizing the excluded region of supersymmetric parameter space. Here we see, for \i > 0, 
the allowed and excluded regions in the S-T plane. The region with T > 0.12 (0.07) is excluded 
in the SUGRA (GM) model. Also, we can read off positive and negative bounds on S. 

The sensitivity of some of the observables is shown in Figs. ^ (a-c) and |TD| (g-1). Here we 
see, for example, 5sin 2 #| ff < -6 x 10~ 4 (-4 x 10" 4 ) and 5A FB (b) > 4 x 10~ 3 (3 x 10~ 3 ) are 
excluded in the SUGRA (GM) model with (i>0. 

The global precision analysis leads to interesting bounds on the superpartner and Higgs 
boson masses. In Fig. [TTJ (a-c) we show the significantly improved limits, > 320, 255, 
230 GeV, in the three models, with \x > 0. In the /i < case these bounds are weakened 
to rriA > 115, 220, and 185 GeV (Figs. [| (d-f)). These limits imply bounds on other Higgs 
boson masses. For example, the region with the lightest possible Higgs boson mass, 62.5 < 
rrih < 71 GeV, is absolutely excluded, in every model and for any sgn(/i)0. The upper bound 
on m/j depends on the cut-off chosen for the squark mass scale. With an 0(1) TeV cut-off 
the upper limit in the GM models, m h < 115 GeV, is stronger by about 10 GeV compared to 
SUGRA due to the relative absence of mixing in the top squark sector. All the Higgs-boson 
mass lower limits primarily result from the B(B — ► X s , -f) constraint on the charged Higgs 
mass, e.g. m H + > 330, 265, 240 GeV, in the fi > case (Figs. |TT| (d-f)). 

The lower bounds on the lightest chargino mass, m~+ > 89, 100, 120 GeV, are seen in 
Figs. [11] (g-i) in the /j > case, together with the corresponding bounds on the lightest 
neutralino mass. The indirect bounds on the squark (specifically ul) and gluino masses can 
be read off of Figs. [U] (j-1) and |8] (g-i), and in the \i > case extend well beyond the direct 
limits. For example, the left-handed up-squark mass is bounded by 285 (355) GeV in the 

15 Interesting bounds on Mq and other supersymmetry parameters have also been obtained by studying 
charge and color breaking minima of the effective potential [Q. However, these bounds are relaxed when one 
considers the possibility that the tunneling time from a false to the true vacuum can exceed the lifetime of the 
universe [p5| . 

16 Table ^reflects a stronger statement of this bound based on more recent LEP results, as discussed in the 
next section. 
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SUGRA (GM) model. 

Both the direct and indirect limits on the third generation squark masses are different from 
the first two generations. The indirect bounds are illustrated in Figs. |9| (a-f) and |12] (a-f). For 
a large range of third generation squark masses, in some cases exceeding 1 TeV, it is possible 
to find excluded parameter space. However, the regions that are absolutely excluded are not 
so significant for these masses. In particular, the light top squark cannot be excluded in the 
SUGRA model (Fig. |T2|(d) ) . However, the heavy (predominantly left-handed) top squark mass 
is significantly bounded by m^ 2 > 275 (355) GeV in the SUGRA (GM) model, with /z > 0. 

The constraints on the slepton masses are shown in Figs. [12] (g-1) and, for the third 
generation, in Figs. | (g~i)- The absolutely excluded bounds for these and other particle 
masses and model parameters are listed in Table [3[ They combine our present knowledge from 
phenomenological constraints^ direct searches and precision experiments. Entries marked by 
an asterisk are determined by the phenomenological constraints or the direct search limits, i.e. 
our global precision analysis does not improve the bound. 

The range of best fit values of m t is confined to the range 168-174 GeV in the parameter 
space with A% 2 < 3.84. Thus, supersymmetry prefers m t in the lower part of the allowed 
Tevatron range. The Z lineshape determination of a s is known to be very sensitive to new 
physics, which in general alters the theoretical prediction for Ri. This happens typically 
through vertex corrections or other corrections to the weak mixing angle. Our analysis reveals, 
however, that in the allowed regions of the investigated models, the fitted a s values are within 
the small range 0.1195 ± 0.0006. Therefore, in the context of minimal gauge- or gravity- 
mediated supersymmetry breaking (including the external constraint 0.118 ± 0.003), 

a s = 0.1195 ± 0.0022 ± 0.0006, 

where the first error is experimental and the second theoretical, i.e. from varying the super- 
symmetry parameters. 

5 Conclusions and Outlook 

Precision data continues to play an important role in discriminating and constraining various 
extensions of the standard model. For example, it is well known that the precision data analysis 
leads to serious constraints on technicolor models. Supersymmetric models can always avoid 
such constraints because the supersymmetric corrections decouple, so the supersymmetric 
models look just like the standard model if the supersymmetric mass scale is large enough 
(~ 0.5 to 1 TeV). As long as the standard model with a light Higgs boson provides a good fit 
to the data (the current fit is excellent), supersymmetric models can as well. 

On the other hand, what if one were to consider a supersymmetric spectrum which includes 
light (O(Mz)) superpartners? Then it is interesting to investigate whether the global fits to 
the data become better or worse as a result of non-negligible supersymmetric corrections. 
This has been the purpose of this paper. To this end, we combined a state of the art standard 
model calculation with the one-loop supersymmetric corrections, and performed global fits to 
the precision observables in three supersymmetric models. We have shown that our global 

17 We refer to electroweak symmetry breaking and Yukawa coupling perturbativity constraints. 
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analysis leads (in many cases) to significantly improved bounds on the parameters and masses 
of the various models. For example, we have shown that for a positive /x-parameter the CP- 
odd Higgs boson mass lower bound increases from ~ 60 GeV to over 320 (230) GeV in the 
supergravity (gauge-mediated) model. For fi < the bound is not as strong, but we are 
able to set a bound, 

m A > 115 GeV, (5) 
for all models under consideration. This has an important consequence: as shown in Figure 21 



of the very recent Ref. p8fl , there is unexcluded parameter space with 60 < nih < 78 GeV 
when 68 < tua < 102 GeV. With our bound (|5|) this part of parameter space is excluded and 
we conclude, 

m h > 78 GeV, (6) 

i.e. the MSSM Higgs bound coincides with the SM one. The rrih bound @ supersedes the 
bound of 71 GeV mentioned in Section |], which has a similar origin, but was based on the 
data presented at the summer conferences [I8~|; for similar bounds with fixed values of tan/? 
see Ref. [|ll[ . 



We set many other interesting bounds, and our main results can be found in Table || and 
Figures |8|-|r2|. We advocate Table |3] as a reference for the current collider and precision bounds 
in an important class of supersymmetry breaking scenarios. 

In summary, significant portions of the parameter spaces of popular supersymmetry break- 
ing scenarios have been shown to be excluded, in particular, for /i > 0. This has immediate 
consequences for the near future of collider experiments. LEP 2 will ultimately be able to dis- 
cover a chargino with a mass up to about 100 GeV. Hence, LEP 2 is very unlikely to discover 
the chargino in the gauge mediation models with fi > 0. The main injector at Fermilab, a 
possible further luminosity upgrade (TEV 33), and the LHC, however, will be able to cover a 
large part of parameter space which is not excluded by either present direct limits or precision 
experiments. 
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Table 2: Results of a global fit to the Standard Model. For each observable, we list the 
experimental result, the best fit result within the SM with the Higgs mass allowed as a free 
parameter, and the pull. 
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Table 3: Limits on model parameters and superpartner and Higgs-boson masses. All masses 
and mass parameters are in GeV (except for A in TeV) and are lower limits. Limits for S, T, 
U, and tan/3 are as indicated (T and U are positive definite). Entries marked by an asterisk 
are consequences of constraints from electroweak symmetry breaking, Yukawa perturbativity, 
or direct searches, and are not significantly improved by our analysis. 
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Figure 4: Supersymmetric contributions to various precision observables in the supergravity 
model. Shown is the effect on each observable (in units of standard deviation) when each 
oblique parameter, S, T, or U, is varied within its possible range in the supergravity model. 
Similarly, we show the full range of the combined effect of all three oblique parameters (labeled 
obi), so that a circle (cross) corresponds to the point with the largest positive (negative) 
oblique shift for that observable. Just below, we indicate the size of the corresponding non- 
oblique corrections at these points. Likewise, the full range of the non-oblique corrections are 
shown (labeled n-obl), and just below the corresponding oblique corrections at the extrema 
of this range are indicated. Finally, we show the range of the full corrections found in our scan 
(labeled full), and at the endpoints of this range we break the corrections down into their 
oblique and non-oblique parts. These are indicated just below, respectively. The non-oblique 
corrections to k are not shown. 
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Figure 5: As Fig. HI for the 5 + 5 gauge-mediated model. 
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Figure 6: As Fig. f|, for the 10 + 10 gauge-mediated model. 
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Figure 7: The range of best fit predictions of precision observables in the supergravity model 
(upper horizontal lines), the 5 + 5 gauge-mediated model (middle lines), the 10 + 10 gauge- 
mediated model (lower lines), and in the Standard Model at its global best fit value (vertical 
lines), in units of standard deviation. 
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Figure 8: Allowed and excluded regions for various observables and supersymmetry param- 
eters in the minimal supergravity model, and the 5 + 5and 10 + 10 gauge-mediated models 
(respectively, from left to right), with /i < 0. The regions which contain allowed points lie 
inside the solid curves, while the regions which contain excluded points are bounded by dashed 
curves. Therefore, the regions inside dahsed curves and outside solid curves are absolutely 
excluded, independent of the values of other parameters. 
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Figure 9: Allowed and excluded /i < regions in third generation squark and slepton mass 
parameter spaces, as described in Fig. Rl Masses are in GeV. 
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Figure 11: Allowed and excluded regions, as in Fig. R with \i > 0. 
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Figure 12: As Fig. [8], in squark and slepton mass parameter spaces, with \i > 0. Masses are 
in GeV. 
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